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ABSTRACT 


The  problem  considered  in  this  paper  is  that  of  finding  the  heat 
losses  from  a  turbulent  gas  flowing  through  a  pipe  due  to  free  and  forced 
convection,  conduction  through  the  solid  insulation  and  external  wall 
radiation.  Furthermore,  along  the  inner  wall,  neither  the  heat  flux  nor 
the  temperature  variation  is  known  in  advance,  but  are  determined  by 
the  losses  from  the  outer  wall.  The  analytical  procedure  used  here  is 
analogous  to  that  first  used  by  Latzko  (1)  and  later  modified  by  Fettis  (7) 
for  a  constant  wall  temperature,  where  eigenfunctions  of  an  approximate 
equation  are  firsi  obtained  and  these  are  then  used  to  obtain  the  solution 
to  the  actual  heat  equation. 

The  numerical  solution  for  the  case  a  =  0.  7  is  presented  in  the  form 
of  graphs  of  axial  temperature  distribution  along  the  pipe  as  well  as  radial 
temperature  profiles  within  the  gas  stream.  The  case  where  the  heat 
transfer  parameter  a  =  0.  7  is  used  corresponds  to  the  specific  dimensions 
and  thermal  properties  of  the  stagnation  chamber  insulation  in  the  ARL 
30-inch  hypersonic  wind  tunnel.  This  tunnel  can  be  briefly  described  as  a 
blow -down  type  facility  which  uses  a  methane -oxygen  fired  zirconia  pebble  bed 
which  heats  the  air  to  the  neighborhood  of  4400°R  prior  to  expansion  to  flow 
Mach  numbers  ranging  from  16  to  22  at  free  stream  Reynolds  numbers 
approaching  10®  per  foot.  The  axial  temperature  distributions  and  radial 
profiles  presented  herein  are  for  a  specific  set  of  conditions  under  which  the 
ARL  30-inch  hypersonic  wind  tunnel  has  actually  operated. 
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NOMENCLATURE 


constant  defined  by  equation  (4) 

constants  determined  from  the  inlet  condition 

parameter  used  in  the  hypergeometric  equation  (27) 

parameter  used  in  the  hypergeometric  equation  (27) 
specific  heat 


parameter  used  in  the  hypergeometric  equation  (27) 
eigenfunctions  for  equation  (21-a) 


“"“SiKS  ZTSSSS*** '» — 

heat  transfer  coefficient  due  to  free  convection 
heat  transfer  coefficient  of  the  gas  in  the  inner  tube 


heat  transfer  coefficient  at  the  inlet  to  the  channel  expressed 
as  unit  thermal  conductance  by  Latzko  *** 


effective  over  all  conduction  coefficient  of  the  pipe  insulation 

SOme  di8tance  away  fro“  ‘he 


thermal  conductivity  of  the  hot  gas  in  the  inner  tube 
heat  transfer  parameter  given  by  equation  (11) 

exponent  in  equation  (4),  depending  on  the  Reynolds  number  range 
Prandtl  number 


parameter  defined  by  equation  (9 -a) 
total  heat  flux 

heat  flux  of  air  tangent  to  the  outer  wall 
heat  flux  at  pipe  inlet 

heat  flux  along  the  radius  of  the  inner  wall  surface 


R 

Re 

r 

r 

2 

r 

i 

? 
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radial  heat  flux  (5lH — ),  positive  in  the  +  r  direction 

hr -ft2 

heat  flux  of  gas  designated  as  region  (2)  in  figure  6 
heat  flux  designated  as  region  (1)  in  figure  6 

heat  flux  due  through  the  solid  layers  from  the  point  tangent 
to  the  outer  layer  of  the  gas 

radius  of  the  inner  tube  surface 


Reynolds  number 
radial  coordinate 

radius  from  center  to  inner  wall  surface  as  indicated  in  figure  6 
radius  from  center  to  outer  wall  surface 


nondimens ional  radius 


subscript  s  refers  to  solid  portion  of  tube,  where  s  and  s  refer 
to  outer  and  inner  solid  surface  layers  respectively 


temperature 

temperature  of  the  gas  immediately  tangent  to  the  outer  wall 

temperature  of  gas  at  the  centerline  of  the  tube. 

temperature  of  the  gas  at  the  inner  wall  designated  as  region  2 
in  figure  6 

temperature  at  the  outer  solid  wall  of  the  cylinder 

temperature  of  the  solid  wall  at  any  point  immediately  adjacent 
to  the  gas  at  the  wall 

reference  temperature  of  the  environment 

effective  thermal  conductance  considered  at  the  radius  of  the 
inner  tube  surface 

velocity  of  the  fluid  stream  in  the  direction  of  flow 
average  velocity 
axia.  coordinate 

nondimens  ional  axial  coordinate 
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o 

heat  parameter  defined  as  — 

7N 

parameter  which  is  proportional  to  the  exponential  decay 
factor  for  the  temperature  in  the  axial  direction 

refers  to  gamma  function 

ratio  of  the  conduction  to  radius  as  defined  by  equation  (7-6) 
under  section  VII 

temperature  difference  between  the  local  temperature  and  a 
reference  temperature 

temperature  difference  at  the  inlet  of  inner  tube 
mean  temperature  difference 

eigenvalues  for  equation  (23) 
kinematic  viscosity 
density 

referring  to  a  summation 
local  shear  stress 
inner  wall  shear  stress 

eigenfunctions  for  the  approximate  differential  equation,  that  is, 
corresponding  to  the  eigenvalues  expressed  by  A 

eigenvalues  for  the  actual  heat  equation  (21 -a) 


I.  INTRODUCTION 


The  problem  of  heat  transfer  in  fully -established  turbulent  flow  in  a 
cylindrical  tube  has  been  investigated  by  many  researchers  (Reference  1 
through  4).  In  all  cases  the  flow  Reynolds  number  is  sufficiently  large  to 
justify  the  assumption  of  negligible  axial  conduction  in  the  fluid.  As  a 
consequence,  the  mathematical  investigations  were  reduced  to  solving  an 
eigenvalue  problem  with  the  wall  temperature  taking  the  form  of  a  step 
function.  Once  this  was  accomplished,  other  boundary  conditions  such  as 
prescribed  heat  flux  or  prescribed  wall  temperature  variation  were  included 
by  the  method  of  superposition,  an  idea  which  was  an  outgrowth  of  Duhamel’s 
theorem.  This  was,  to  the  author’ s  knowledge,  first  pointed  out  by  Tribus. 
However,  there  are  certain  situations  in  which  the  use  of  the  superposition 
integral  may  not  be  practical.  Such  is  the  case  in  a  cross-flow  heat  exchanger 
when  the  wall  heat  transfer  rate  is  related  to  an  external  heat  transfer 
coefficient.  In  this  case,  neither  the  heat  flux  nor  the  temperature  variation 
is  known  in  advance  and  the  flux  must  be  matched  with  the  wall  temperature 
variation.  To  accomplish  this  one  needs  to  know  the  eigensolutions  and  the 
coefficients  quite  accurately  near  the  entrance. 

The  practical  engineering  use  for  the  steady  state  heat  transfer  analysis 
and  numerical  example  presented  here  is  that  of  determining  the  heat  losses 
from  a  turbulent  gas  flowing  through  an  insulated  pipe.  A  cross-section  of  the 
pipe  considered  in  this  study  is  shown  in  Figure  1.  Consideration  is  given  to 
finding  heat  losses  due  to  all  three  modes  of  heat  transfer  with  a  further 
stipulation  that  the  surface  temperature  of  the  chamber  is  not  held  constant. 

In  the  first  portion  of  the  paper,  consideration  is  given  to  solving  the 
general  heat  transfer  problem  for  turbulent  fluid  flow  through  cylindrical 

tubes.  The  analytical  approach  for  solving  the  problem  is  given  in  the 

Manuscript  released  by  the  author  April  1965  for  publication  as  an 
ARL  Technical  Report. 
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following  three  sections.  The  remainder  of  the  paper  gives  a  numerical 
example  for  the  specific  case  corresponding  to  the  various  thermal 
conditions  in  the  stagnation  chamber  of  the  ARL  30-inch  hypersonic  wind 
tunnel. 
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H.  STATEMENT  OF  THE  PROBLEM 


The  problem  is  stated  as  follows:  given  the  initial  temperature 
distribution  (with  respect  to  the  radial  distance)  of  a  fluid  entering  a 
cylindrical  tube  through  which  heat  transfer  occurs  by  means  of  a  fixed 
wall  conductance,  the  ensuing  temperature  distribution  is  to  be 
determined. 

Since  the  accent  of  the  investigation  is  on  this  particular  boundary 

condition,  the  conventional  assumptions  of  constant -property  (including 

density)  fluid  are  adopted.  The  flow  Reynolds  number  is  in  the  order  of 
5 

10  and  the  length -diameter  ratio  is  large  so  that  very  little  error  is 
introduced  by  assuming  the  existance  of  fully -established  flow  beginning 
at  the  entrance  section.  Ignoring  axial  conduction,  the  steady-state 
energy  equation  with  circular  symmetry  is 

3T  i  a 

-p  CP  U  ax  =  7  7r  (rqr)  (1) 

where  qr  is  the  radial  heat  flux,  positive  in  the  +  r  direction.  On  the 
inner  wall  r  =  R,  the  flux  is  equal  to  the  thermal  conductance  UR 
multiplied  by  the  temperature  difference  between  the  inner  wall  and  a 
reference  temperature,  Trej  .  Thus,  at  r  =  R: 

«R  =  UR  <T-Tref>  •  <2) 

Uj^  includes  contributions  from  the  wall -thickness,  convection 
and  external  wall  radiation  effects.  It  is  assumed  that  is  effectively 
independent  of  axial  distance  x.  At  the  entrance,  the  fluid  temperature 
has  its  initial  value  9  (r),  that  is  at  x  =  o 

8  =  8q  (r) .  (3) 

The  equations  (1),  (2),  and  (3)  for  "q"  constitute  the  present  problem. 


IE.  ANALYTICAL  METHODS 


Currently  there  are  four  different  approaches  toward  the  solution 
of  the  stated  problem.  They  can  be  classified  into  two  main  categories; 

(i)  eigenvalue  approach  represented  first  by  the  work  of  Latzko  (1)  and 
later  modified  by  Fettis  (7),  Sleicher  and  Tribus  (2)  and  Beckers  (3), and 

(ii)  boundary  layer  approach  typified  by  the  work  of  Deisler  (4)  and 
Sparrow  (5).  As  a  general  conclusion,  the  analysis  of  (4)  and  (5)  based 
on  the  integral  method  of  the  boundary  layer  equations  gives  a  gross 
picture  without  a  high  degree  of  precision  in  its  numerical  results. 
Among  the  eigenvalue  approach  techniques,  there  are  three  different 
ways  of  handling  the  turbulent  velocity  distribution  and  the  eddy 
diffusivity  distribution.  Latzko  employed  a  simplified  1/7  power  law 
closed  form  distribution  for  the  velocity  and  the  total  thermal  diffusivity 
(turbulent  plus  molecular).  Sleicher  and  Tribus  used  experimentally 
measured  velocity  distributions  and  eddy  diffusivities.  This  resulting 
procedure  was  to  solve  the  problem  for  each  Reynolds  number  and 
Prandtl  number. 

To  use  the  method  given  by  Beckers,  one  needs  to  divide  the  flow 
into  three  different  regimes;  viscous  sublayer,  buffer  layer  and 
turbulent  core.  In  fact,  he  solved  three  sets  of  differential  equations 
with  three  solutions  which  are  joined  together  with  temperature 
continuity  only. 

The  approach  of  Sleicher  and  Tribus  appears  to  be  the  best 
procedure,  but  the  numerical  work  for  a  general  treatment  is  also  the 
most  time-consuming.  The  second  best  is  the  Latzko  method  which 
can  be  modified  for  a  different  Reynolds  number  range  by  using 
different  power -law  exponents  in  his  original  analysis.  In  view  of  the 
flexibility  afforded  by  the  Latzko’ s  formulation,  a  portion  of  the  latter 
approach  is  adopted  in  the  present  problem. 
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IV.  ANALYSIS 


Latzko's  analysis  is  based  on  the  premise  of  the  power -law  distribution 
of  velocity  near  the  wall  which  works  quite  well  in  that  vicinity.  However, 
since  the  velocity  profile  becomes  relatively  flat  near  the  tube  axis,  he 
expressed  the  velocity  distribution  with  the  following  expression, 


u  = 


n 


i 


Equation  (4)  reduces  to  u  ~  6  n  when  r  is  close  to  R,  that  is  6  =  R  .  r  «  o. 
The  coefficient  "A”  is  related  to  the  mean  flow  by  a  constant  factor  dependent 
on  n. 


For  fully  established  flow,  the  total  (molecular  and  eddy)  shear  is  a 
linear  function  of  the  radius.  If  equation  (4)  is  differentiated,  one  obtains  the 
following, 


In  general  n  can  take  on  values  of  6,  7,  8  and  9  depending  on  the  various 

Reynolds  number  ranges  considered;  hence  the  velocity  derivative  exhibits  a 

singularity  at  r  =  R,  indicative  of  the  fact  that  the  power  law  expression 

cannot  be  used  at  the  wall.  However,  the  quantity  on  the  left  hand  side  of 

equation  (5)  remains  regular  at  r  =  R  which  must  be  equal  to  the  shear  ratio 

(  t/t  ),  where  t  is  the  wall  shear.  Accordingly,  Latzko  expressed  the 
R  R 

shear  ratio  as, 
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t  _  -nR 

T  A 
R 


8u 

dr 


Further,  he  used  the  Reynolds  analogy  for  a  Prandtl  number  equal  to  one  and 
obtained  an  expression  for  the  heat  flux  using  n  =  7, 


In  this  manner  the  temperature  derivative  is  permitted  a  singularity  at  r  =  R, 
but  the  heat  flux  as  defined  by  (6)  is  regular  at  r  =  R.  Quantitatively,  Latzko 
used  the  following  expression  for  heat  flux, 


3/4 

-0. 199  V 

3/ 18 

(2R) 


1/4 

v  pCp 


d0 

dr 


(6-a) 


and  in  the  presently  used  notation,  qr  is  expressed  as 


With  the  constant  A  in 
expressed  as 

8  „ 
u  =  -  V 

7 


equation  (4)  evaluated  for  n  = 


7,  the  velocity  is 


Introducing  the  non-dimensional  variables  r  =  l-i  x 

,  x  Ir / 

temperature  difference  as  0  ,  equation  (1)  reduces  to 


I  -j  and  the 

Irj 


6 


(9) 


p(l-r2)*^  =  i  L  \r  (l-f2)’/7  M 

0x  f  3f  L  3r 

where  p  =  f—  2  — *  |  Re  ‘/4 
1 7  x  0. 199  J 

The  boundary  condition  at  r  =  R  is  q^  =  UR(T-Tref>  or 


(9 -a) 


'  °'  199  k  P  R 

»  < 

i  -  f2 

* 

e/7  1 

2  1  i7  g  r  e 

*  < 

R 

for  r  =  1, 


80 

dr 


(T  -  Tref) 


(10) 

(ID 


The  temperature  distribution  at  x  =  o.  becomes 

0  =  3  Q  (r)  (12) 

where  6  represents  the  temperature  difference  at  the  inlet  of  the  pipe. 

A  second  condition  must  be  satisfied  which  demands  that  the  centerline 
temperature  be  finite.  This  second  condition  will  be  given  after  another 
change  of  variables  is  made  and  substituted  into  equation  (9). 


A  solution  to  equation  (9)  is  assumed  to  be  of  the  form, 
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Substitution  of  equation  (13)  into  equation  (9)  yields  the  following 


2  /1  1/7  ^  Id 

■p  (1  -  r)  /  g  =  -  — 

r  dr 


-2\8 


r  (1  -  fa) 


/7  dgl 

dr  J 


Introducing  the  change  of  variables 


z7  =  1  -  r 


-  u 


equation  (14)  becomes 


d_ 

dz 


(1  -  27  ) 

dz 


+ 


g  =  o 


and  equation  (10)  is  written  as: 


( 16-a) 


+  g  =  o  (16-b) 

z  =  o 

The  second  condition,  mentioned  previously,  is  to  satisfy  equation  (16-a)  at 
z  =  1  or  the  centerline.  Because  of  the  singularity  at  z  =  1,  it  is  necessary 
to  require  that  g  be  finite  at  z  =  1.  The  system  described  by  (16-a)  with 
proper  boundary  conditions  is  a  regular  Sturm-Liouviile  system  solved  as 
an  eigenvalue  problem.  The  o  thogonality  relation  given  below  is  readily 
established, 


.  2_  4 

7N  dz 


A  z?  gn  ^  gm  ^  ^  =  0  ior  m  ^  n 
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The  initial  conditions  are  satisfied  by  putting 


e=lAn  exp |  xj  g„  (z)  , 

where  the  coefficients  A„  are  determined  from 

n 

*  =  Jo  gp  gn  (z)  d Z 
n  ri 

J  o  z7  gn  (z)  d z 

and  6q  is  defined  at  x  =  o,  or  60  =  6  (z,  o). 


Let  = 


49 

4 


(20) 


then  equations  (16-a)  and  (16-b)  respectively  reduce  to 


+  U)  z  g  =  o 


(21-a) 


Evaluating  at  z  =  o,  and  for  convenience  in  later  computing  the  equation 


involving  gamma  functions,  letting  —  =  a ,  we  have 

7N 


( 


dg1 

dz 


+  g 


z  =  o 


=  o 


(21-b) 


z  =  o 


Solutions  of  (21-a),  which  satisfy  the  condition  stipulated  at  z  =  1  and  also  of 
equation  (21-b)  will  c^ist  only  for  discrete  values  of  the  parameter  as 
indicated  in  equation  (21-a).  Once  these  eigenvalues  have  been  found,  the 
complete  solution  of  equation  (9)  will  then  be  of  the  form 
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(22) 


ZAnKn  (z)exp 


Before  considering  the  solution  of  (21-a),  we  wish  to  consider  the 
solution  of  an  auxiliary  equation  which  is  similar  in  appearance,  although 
not  identical,  to  equation  (21-a),  that  is: 


—  (1  -  z7  )  =  -  X  z“  g 

dz  dz 


} 


(23) 


A  change  of  variable  z7  =t  converts  equation  (23)  to  the  following 


t  (1  -  t)  +  [  -  -  --  \  +  —  g  =  o 

dta  17  7  J  dt  49 


(24) 


A  similar  equation  results  if  we  consider  still  another  change  of  variables, 
this  being 

y  =  1-t 
and  dy  =  -dt 


(25) 


Then  equation  (24)  becomes 


(1  -y)  +  (l  -I3-y\  &  +  ±g  =  o 

dy*  \  7  J  dy  49 


(26) 


Equations  (24)and  (26)  are  in  the  form  of  the  hypergeometric  equation 


y  (i  -  y) 


c  -  (a  +  b  +  1)  y 


^  -  a  b  g 
dy 


o 


(27) 


Equation  (24)  is  useful  when  obtaining  solutions  for  a  =  o  and  ®.  When 
a  =  o,  then  9=o  which  corresponds  to  Latzko*  s  solution  to  the  heat  transfer 
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problem,  also,  when  q  =  ®,  the  solution  corresponds  to  the  case  =  o. 

For  the  case  of  q^  =  o,  a  meaningful  expression  for  the  temperature  at 
x  =  o  must  be  prescribed,  that  is,  0Q  cannot  be  independent  of  the  radius. 
Also,  a  solution  is  sought  such  that  the  initial  inlet  expression  for  0Q  is 
successively  "damped  out.  ”  For  all  other  values  of  the  parameter  "a", 
we  must  obtain  solutions  to  equation  (26). 

With  the  aid  of  equation  (27),  comparison  of  coefficients  a,  b  and  c  can 
be  made  which  are  used  in  equations  (24)  and  (26).  For  equation  (24)  we  obtain 
the  following; 


7 


7 

while  for  equation  (26),  we  obtain 


(28) 


7 

c  =  1 

For  the  general  case,  equation  (27),  which  includes  a  =  o  and  a 
linearly  independent  solutions  are 


(29) 


®,  two 
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(30-a) 


gl  =  F  (a,  b,  1,  y) 

gt  =  Fj  (a,  b,  1,  y)  1  n  y  (30-b) 

but  because  of  the  requirement  that  g  be  finite  at  z  =  1,  only  the  solution 
(30-a)  is  admissible.  Hence,  except  for  a  multiplicative  factor,  the  desired 
solution  of  (27)  is 


g  =  F  (a,  ^  -  a,  1,  y)  (31) 

The  value  of  "a",  and  therefore  A ,  is  still  undetermined,  and  must  be  so 
chosen  that  the  boundary  conditions  are  satisfied. 

From  the  known  formula  for  the  derivative  of  the  hypergeometric  function 
called  "  F” 


—  =  p(a+  1,  b+  1,  c+  1,  y) 
dy  c 


the  derivative  of  (31)  can  thus  be  written  as 


(32) 


dg  = 

dy 


a)  F  (a  +  ~1 


-  a,  2,  y 


(33) 


At  y  =  1,  5S  is  singular.  However,  the  boundary  condition  is  given  in  terms 
d^  dy 

of  and  since  y  =  1  -  z  ,  we  find  that 
dz 


dg  =  dg  djr  =  .  7  z«  dg  (34) 

dz  dy  dz  dy 
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Equation  (33)  can  now  be  written  as 


dy 


[-  - a1 

p 

l7  i 

\ 

13 

7 


(35) 


Now  the  appropriate  general  expression  for  the  function  F  is  the  following: 


F  (a,  b,  c,  y)  =  T-^  r-(c-  — a-_  b)  F  (a,  b,  a  +  b  -  c  +  1,  1  -  y) 

r  (c  -  a)  r  (c  -  b) 


+  (1  -  y)C  a  k  —  ^  F  (c  -  a,  c  -  b,  c  -  a  -  b  +  1,  1  -  y) 

T  (a)  T  (b) 

(36) 


where  y  =  1  -  z7  . 
Hence, 

13 


F  (a  +  1,  ~  -  a,  2,  y) 
7 


T  (2)  T  (-y) 

r  (1  -  a)  r  (-  +a) 
7 


F  (a  +  1,  --a,  1^,  z7) 


+  z 


-6 


r  (2)  r 

A 

7 

r  (a  +i)  r 

PH 

F  (1  -  a,  i  +  a,  -  ,  z7  ) 
7  7 


(37) 


Substituting  (37)  into  (35),  one  obtains  the  following: 

r  (2)  r  (-■§-) 


F  (a  +  1,  H 


(1  -a)r  (i  +  a) 
7 


-a;  z’) 


+  z 


r  (2)  r  (■§■) 


r  (a  +  l)  r  (— -  a) 

7 


F(1  -  a,  i  +  a,  i,  z7  ) 


(38) 
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But 


dg 

dy 


dg=  dg  djr  _  _7z« 

d z  dy  dz 

hence, 


lim  &  =  -  7a  ( *  -  a) 
z~+o  dz  7 


r  (2)  r  (-fj 
r  (a  +  1)  r  (—  -  a) 

7 


F  (1  -  a,  i  +  a,  i,  o) 

7  7 

(39) 


At  the  origin, 

F(l  -  a,  -  +  a,  0)  =  l 
7  7 


Therefore 

lim  dg=  -7  r  <2>  r  4 

z— o  dz  r  (a)  r  (-  -  a) 

7 

At  z  =  o,  equation  (31)  becomes 

g  (o)  =  F  (a,  -  -  a,  1,  1) 

7 


(40) 


where 

F  (a,  b,  c,  1)=  r  i£>  r  (c  -  a  -  b)  (42) 

T  (c  -  a)  T  (c  -  b) 


We  now  obtain  from  equation  (21 -b),  with  the  aid  of  equations  (41)  and 
(42),  the  relation 


7a  r  ($) 
r  (a)  r|(|  -a)j 


r 


r  (1  -  a)  r  (i  +  a) 


(43) 
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Thus  admissible  solutions  occur  when  "a"  is  such  that 


r  (-i-)  r  (a)  r  a) 

a  =  - 1 - 1 -  (44) 

7  r  (i  -  a)  r  (i  +  a)  r  (p 

The  values  for  A  corresponding  to  values  of  "a”  which  satisfy  equation 
(44)  can  then  be  obtained  from  the  formula 

X  -  -  49  ab  ==  -  49a  (-  -  a)  (45) 

7 

For  values  of  a  other  than  zero  or  infinity,  the  eigenfunctions  of 
equation  (23)  may  be  obtained  in  the  following  form  from  use  of  equation  (31): 

«,W=  An  F  (a„ ,  8  -an,  1,  l-z7)=  An  ^  (z)  (46) 

The  functions  f  (x)  are  orthogonal  with  respect  to  z5  as  weight  factor, 
that  is 

/o  z*  (Z)  fn  (z)  dz  =  o,  if  m  /  n  (47) 

and 

So  z’  fn*  1,2  =  So  z>[F(an’°  ~  an  ’  ^  1  >  1  dz  =  Cn  <48> 

For  application  to  equation  (21 -a),  it  is  convenient  to  define  a  new  set  of 
functions  $n  (z),  for  the  solutions  of  the  approximate  equation  (23)  such  that 

L  (z) 

®n  (z)  =  “  -  (49) 

n  c.vrr 
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and 


A  „  /*  z5  9  (z)  *  (z)  dz  =  6 

n  •'o  m  n  '  ' 


mn 


(50) 


where 


6 _ =  0,  if  m  i  n;  and 

mn  ’  '  ’ 


6 _ =  1,  if  m  =  n 

mn 


Now,  if  we  assume  a  solution  of  the  form 
S 

g  (z)  =  z  Bn  *n  ^ 
n=  1 


(51) 


and  substitute  this  expression  into  equation  (21-a),  we  obtain 


s  r 


(1  -  z 
dz 


d  *  1 

7  )  - HL  +  u>  z7  4> 

dz  n  n 


=  0 


and  because  the  $n  's  satisfy  equation  (23),  this  becomes 


5  f  1 

2  Bn  Xn  z’  *n  (z)  *  wn  z’  *n  (z) 

n=1  L  J 


=  0 


(52 


(53) 


Thus  equation  (51)  is  the  desired  solution  to  equation  (21-a). 

To  find  the  coefficients  Bn ,  we  require  that  the  left  hand  side  of  this 
equation  be  orthogonal  to  the  4>n  for  n  =  1,  . . . ,  S  . 

This  implies  the  following, 

S 

Z  Bn  ( 6  -  A  )  =0;  m  =  1,  . . . ,  S  (54) 

u  n  mn  n  mn 

n=  1 
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where 


A 


mn 


/*  z7  $  (z)  $  (z)  dz 

Jo  m  7  n  v  7 


(55) 


The  characteristic  equation  of  the  system  (54)  is 


6 

mn 


0 


(56) 


and  the  roots  of  equation  (56)  give  approximations  to  the  first  S  eigenvalues  of 
equation  (21-a).  Further,  with  each  u>n  there  is  associated  a  set  of  B's,  and 


S 

linear  combination  Y  ®  (z)  is  an  approximation  to  the  corresponding 

n=  i 

eigenfunctions  of  equation  (21-a). 
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V.  NUMERICAL  SOLUTION  FOR  THE  SPECIAL  CASE  a  =  0.  7 


The  specific  value  of  a  to  be  introduced  at  this  point  is  taken  to  be 
that  corresponding  to  the  conditions  in  the  airstream  at  the  inlet  to  the 
stagnation  chamber  of  the  ARL  30-inch  hypersonic  wind  tunnel.  For  the 
special  case  treated  in  this  report,  a  specific  set  of  conditions  has  been 
chosen  which  corresponds  to  tunnel  operation  at  approximately  Mach  20. 
Therefore,  the  conditions  to  be  used  in  the  subsequent  numerical  analysis 
to  describe  the  gas  at  the  pipe  inlet  are  given  below,  while  the  derivation 
of  the  coefficients  for  free  convection,  outer  surface  radiation  and 
conduction  from  insulation  will  be  shown  in  the  next  section.  The  conditions 
at  the  pipe  inlet  are  the  following: 


^inlet  =  Pressure  a*  PiPe  in*e*  =  2000  psia 


T, 


inlet 


m 

A 

P 


=  temperature  at  pipe  inlet  =  4000°R 
=  air  mass  flow  =  0.  53  lb  mass/ sec. 

-3  2 

=  inner  cross-sectional  area  of  pipe  =  3.  068  x  10  ft. 

=  dynamic  viscosity  =  480  x  10"  *  lb  mass 


(sec)  (ft) 


=  specific  heat  =  0.  30  B  Tu 


(lb)  (°R) 

p  =  air  density  =  1.  35  lb/ft^ 

V  =  average  airstream  velocity  =  128.  ft/ sec 

Re  =  Reynolds  number  =  2.  2515  x  10^ 

For  a  =  0.  7,  the  first  four  values  of  "a"  which  satisfy  equation  (44)  are 


a  =  1. 2862 

1 

a  =  2. 4374 

2 

a  =  3. 5195 

3 

a  =  4. 5695 

4 

Using  these  values  of  "a",  one  can  now  obtain  from  equation  (45)  the 


(57) 


values  of  An  corresponding  to  the  values  a^  ---a 


n 


Also  the  eig  enfunctions 
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f  for  equation  (23)  may  be  obtained  by  making  use  of  equation  (46).  Making 
use  of  the  orthogonality  relationship  given  by  equation  (48),  we  can  obtain 
the  coefficients  Cn  which  are  used  in  equation  (49)  to  obtain  the  *nf  s. 

We  can  now  proceed  to  obtain  the  eigenfunctions  gn  (z)  by  use  of  equation  (51). 
But  before  equation  (51)  can  be  solved  explicitly  for  gn  (z)  the  eigenvalues 
ujn  and  coefficients  Bn  must  be  found.  These  can  be  obtained  by  solving 
equation  (53).  Because  of  the  orthogonality  requirement  stipulated  by 
equation  (53),  use  of  equations  (54)  and  (56)  must  be  made  to  obtain  the 
coefficients  Br  and  u;n  .  The  roots  of  equation  (56)  give  approximations 
to  the  first  MS"  eigenvalues.  Making  use  of  these  eigenvalues  and 
equation  (54),  we  can  obtain  the  coefficients  Bn>  With  the  coefficients  Bn 
and  $>n  calculated,  we  can  now  find  the  eigenfunctions  gn  from  equation  (51). 
Before  the  complete  solution  for  6  can  be  expressed,  the  constant  ^  must 
be  solved  explicitly  by  making  use  of  equation  (19).  Having  obtained  the 
coefficients  Ar,  we  can  now  write  the  complete  solution  for  9  as  given  by 
equation  (18). 

Using  the  values  of  "a",  and  the  procedure  as  outlined  in  the  above 
paragraph,  we  can  now  obtain  the  eigenfunctions  gn  (z)  of  the  original 
equation  (21  -  a)  in  terms  of  the  eigenfunctions  4>n  (z)  of  equation  (23). 

Making  use  of  the  IBM  1620  and  7090  high  speed  computers,  there  resulted 
the  following  cun  and  gn  (z)  when  a  =  0.  7. 

w  =  30. 689 

g  =  0. 99830  $  +  0. 05 . 50  4>  -  0.  00988  9>  +  0.  00418  * 

1  12  3  4 

<*>t  =  232. 436 

g  =  -0.05431  «  +  0.97860  *  +  0.19690  *  -0.02516  * 

1  13  3  4 

u)  =  581. 542 

3 
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r  =  0.01776  *  -  0.17540  4  +  0.92100  4  +  0.34750  * 

63  1234 

u)  =  1196.  131 

4 

g  =  -0.01249  ft  +  0.09104  ft  -  0.33610  ft  +  0.93730  ft  (58) 

64  1  2  3  4 

In  the  special  case  whei  a  6 Q  ,  the  initial  temperature  distribution,  is 
taken  as  a  constant,  the  solution  expressed  analytically  by  equation  (18)  can 
now  be  written  as  the  following: 


1.  286  exp[l£ 
j)R 


rx)[ 


0.99830  4  +  0.  05750  *  -  0.00988  *  +  0.00418  4 

i  a  3 


1  = 
*0 


•] 

-0. 399  exp|-Pi  xU.q.  05431  4j  +  0.97860  *a  +  0.19690  43  -  0.02516  44j 

} 

•} 


+  0.  265  exp 


-0.  217  expfl^ 


(?*L  ]  1° 

tj  [-° 


.01776  4  -  0.17540  4  +  0.92100  4  +  0.34750  4 

1  2  3 


.01249  4  +  0.09104  4  -  0.336100  4  +  0.93730 

’  1  2  3 


where  the  n 
pR 


P*  are  as  follows: 


P?  =  0. 3906 

pR 

IL  =  2.9583  (  60) 

pR 

'll  =  7. 4015 

pR 

PJ  =  15.2236 

pR 

The  eigenfunctions  ftn  are  shown  as  a  plot  of  z  versus  ft  in  Figures  2,  3,  4 
and  5. 
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VI.  CALCULATION  OF  TEMPERATURE  DISTRIBUTIONS  AS  FUNCTIONS 
OF  RADIUS  AND  AXIAL  DISTANCE  DOWNSTREAM 

1.  Introduction 

Once  the  explicit  solution  for  temperature  difference  ”d  "  is  obtained 
from  equation  (59),  we  can  proceed  to  find  actual  temperatures  as  functions 
of  radius  and  distance,  ft  is  more  useful  to  the  engineer  to  have  the  results 
expressed  as  temperatures  rather  than  temperature  differences. 

To  obtain  the  temperature  as  a  function  of  location,  all  three  modes 
of  heat  transfer  should  be  considered  simultaneously.  In  view  of  the 
different  laws  governing  the  three  modes,  namely:  conduction,  radiation, 
and  convection,  an  exact  solution  is  extremely  difficult  if  not  impossible. 

In  the  situation  of  heat  transfer  through  pipe  insulation,  one  can  obtain 
reasonable  results  by  treating  all  the  above  named  processes  independently. 
We  wish  to  give  a  simplified  engineering  approach  for  the  solution  of  such 
a  problem  by  making  use  of  the  data  obtained  in  the  previous  section.  The 
stipulation  that  the  temperatures  and  temperature  gradients  be  continuous 
between  any  two  surfaces  exchanging  heat  is  made  in  obtaining  the  final 
results. 

2.  Simplified  Engineering  Analysis 

By  requiring  that  the  temperature  and  temperature  gradient  be  continuous 
between  any  two  surfaces  exchanging  heat,  we  are  able  to  obtain  analytical 
expressions  for  the  centerline,  inner  and  outer  wall  temperature 
distributions  explicitly.  The  geometry  of  the  cylindrical  chamber  cross 
section  being  considered  herein  is  made  up  of  n  -  concentric  layers  of 
different  insulating  materials  of  various  thicknesses  as  shown  in  Figure  6. 

The  procedure  for  matching  the  regional  temperature  gradients 
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and  temperatures  between  any  two  surfaces  exchanging  heat  will  be 
described  in  the  next  sub-section. 

3.  Matching  Techniques 

With  a  known  expression  for  ”0  ”,  and  a  given  reference  temperature, 
the  appropriate  formulae  needed  for  obtaining  surface  temperatures  will 
now  be  derived.  At  the  inner  wall  shown  in  figure  six  as  region  (2),  we 
proceed  to  equate  the  temperatures  and  temperature  gradients  of  the 
fluid  tangent  to  the  solid  wall  and  the  wall  itself.  This  can  be  expressed 
analytically  by  the  following; 


AV  =  A%2 


BT  =  BT 

gW _ 


(6-1) 


Aqgw  +  BTgw  =  A«s,  +  BT* 

Now,  the  heat  flux  q  can  be  expressed  by  the  following  equation, 

gw 

Vv  '  hgR  (Tc  "  V  ■  hg  (Tc  ■  T2^  ■  hgR0 
with  h  ^  defined  in  appendix  B. 

Rewriting  equation  (6-2)  as, 


(6-2) 


v  -  hg  (Tc  ’  V  =  0  <6’3> 

and  also  after  subtracting  BT  from  both  sides  of  equation  (6-1),  one  obtains, 

AV +  B  (Tew  *  V  =  a%2  +  B  Tc)  (6'4) 

Equating  coefficients  A  and  B  from  equation  (6-4)  to  those  in  equation 
(6-3),  one  finds  that, 

A  =  1  and  B  =  -h 

o 
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.  +  h  (T  -  T  )  =  o 

lSa  g  c  2 


Therefore, 
li 

From  the  solution  of  the  steady  state  conduction  problem,  we  have  the 
following  familiar  flux  equation, 


(6-5) 


zJL 

r 


(T  -  T  ) 
2  1 


(6-6) 


From  equations  (6-5)  and  6-6),  we  obtain  the  following  important  relationship 
between  forced  convection  and  overall  conduction, 


1 1 


(T  -  T  )  =  h  (T  -  T  ); 
a  i  g  C  a  ’ 


or  at  region  (2)  we  have  the  temperature  difference  "0  "  expressed  as 


(6-7) 


(T.  -T)=  (T  -T). 


a  i  '  t j  '  c  a 

The  same  analogy  can  be  used  to  obtain  expressions  for  temperature 
relationships  between  heat  conduction,  free  convection  and  radiation  at  the 
outer  surface  of  the  pipe.  Equating  the  gradients  and  temperatures  of  the 
outer  surface  of  the  pipe  with  a  given  reference  temperature  some  distance 
away  from  the  surface,  one  can  use  the  equations  for  continuous  flux  and 
temperature  similar  to  equation  (6-1) 

A<ls,  = 


(6-8) 


BT  =  BT 


(6-9) 

i  a 

The  flux  q  immediately  tangent  to  the  outer  surface  of  the  cylinder,  can  be 


expressed  as, 


qa  *  <Ta  -  Tre{)  =  o 


(6-10) 
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Performing  the  same  type  of  manipulation  here  as  for  obtaining  equation 
(6-5),  we  obtain  the  formula  for  external  wall  flux, 

fls.  -  h  *  Tref>  =  0 


But  we  know  that  the  expression  for  qgi  is 


lSl 


keff 


(T  -  T  )  = 
2  1 


3) 

r 


(T  -  T  ) 
2  1 


(6-12) 


Substituting  for  q  from  equation  (6-12)  into  (6-11),  we  obtain  the  following 

S 1 

relationship, 

<W=  V”  (T.-Tref)-  (6‘13) 

If  equations  (6-8)  and  (6-13)  are  now  combined,  we  obtain 


r  h  (T  -  T  )  =  r  h  (T 
2  g  '  c  r  i  ' 


ref 


(6-14) 


or 


r  fi 


(6-15) 


We  wish  to  express  the  temperatures  T  ,  T  ,  and  T  as  function  of  the 

X  2  C 

parameters  h^,  h  and  r\  as  well  as  the  radii  rj ,  r a,  and  temperature 
difference  6 .  These  relationships  in  algebraic  form  are, 


(6-16) 

(6-17) 


t  =  e  +  t 

C  a 


(6-18) 


24 


After  some  further  manipulation  of  equations  (6-16),  (6-17)  and  (6-18), 
we  obtain  the  following  useful  formulae  for  the  temperature. 


(6-19) 


(6-20) 


and  denoting  the  portion  in  the  brackets  as  coefficient  "c  "  we  have 


T  =  T  ,  +  c  0 
c  ref  i 


(6-21) 


We  can  also  express  the  temperatures  and  T  as  follows, 


T 

2 


(6-22) 


where  c  is, 
2 


) 


and  for  T  the  following : 

T  =  T  -  c  9  (6-23) 

i  ref  3 

with  c3  expressed  as, 


To  obtain  actual  temperatures,  the  coefficients  h,  ke^,  etc.  must  be  defined 
and  expressed  analytically.  For  the  special  case  where  the  value  of  the  heat 
transfer  parameter  a  was  0.  7,  the  formulae  for  the  coefficients  mentioned 
above  are  given  in  the  Appendix. The  results  obtained  for  the  case  a  =  0.  7 
are  presented  (Figs.  9-11)  as  plots  of  axial  temperature  distribution  along 
the  pipe.  Some  radial  temperature  proiiles  at  specific  axial  locations  are  also 
presented  ( Figs.  12-15). 
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VII.  CONCLUDING  REMARKS 


1.  The  plot  of  An,  u>n  versus  q,  shown  in  Figure  8,  indicates  that 
the  exponential  decay  of  the  temperature  along  the  channel  increases 
steadily,  at  first  rather  slowly  up  to  q  =  0.  5  then  quite  rapidly  until 
the  asymptotic  value  (a  =  ao)  is  reached.  This  trend  is  the  same  for 
equations  (21-a)  and  (23). 

2.  The  accuracy  in  approximating  the  uniform  radial  temperature 
distribution  at  the  pipe  entrance  is  dependent  upon  the  number  of 
eigenfunctions  used.  With  only  four  eigenfunctions,  the  accuracy  was 
poor  near  the  wall  but  excellent  near  the  centerline.  However,  the 
accuracy  of  the  radial  distribution  at  the  wall  improves  considerably 
at  subsequent  axial  positions. 

3.  For  the  numerical  example  chosen,  the  fully  developed  thermal  field 
is  attained  when  the  gas  has  passed  through  136  tube  diameters.  Also, 
the  centerline  temperature  dropped  rapidly  from  its  maximum  value  in 
96  tube  diameters,  then  asymptotically  approached  a  minimum  value, 
while  the  corresponding  radial  distribution  became  uniform. 

4.  In  the  axia1  direction  of  the  fluid  downstream,  the  temperature 
distribution  becomes  less  and  less  dependent  upon  the  higher  eigenfunctions. 
This  implies  that  to  obtain  a  first  approximation  to  a  very  complicated 

heat  transfer  problem,  one  need  only  to  make  computations  using  the  first , 
or  at  most,  first  and  second  eigenfunctions. 

5.  In  using  the  eigenvalue  -  eigenfunction  approach  when  obtaining 
approximations  to  the  solution  for  the  heat  transfer  equation,  some  inherent 
error  exists  in  the  (n  -  1)  and  n^  terms  of  the  hypergeometric  series.  The 
effect  of  the  error  from  approximating  the  last  two  terms  is  noticed  only 
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when  approximating  a  uniform  temperature  function  at  the  entrance 
section. 

6.  The  axial  and  radial  temperature  profiles,  as  presented  in 
Figures  9  through  15  inclusive,  may  be  considered  reasonably  accurate 
for  engineering  purposes.  One  reason  for  any  inaccuracy  in  applying 
these  results  would  be  that  near  the  entrance  section  the  velocity  profile 
may  not  be  fully  established  at  the  actual  thermal  t  .trance  section. 

The  analysis  contained  the  assumption  that  the  velocity  profile  was  fully 
established  at  the  thermal  entrance  section. 
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APPENDIX 


In  the  heat  transfer  from  the  outer  wall  of  the  cylinder  to  a 
surrounding  gaseous  medium,  e.  g.  air,  free  convection  and  external 
radiation  are  the  contributing  factors  in  calculating  the  temperature  at 
the  outer  surface  of  the  cylindrical  wall.  The  two  contributions  to  the 
heat  transfer  due  to  free  convection  denoted  as  hc  and  that  due  to 
radiation  defined  as  hr  are  additive  as  a  first  approximation.  Therefore, 
the  heat  flux  q  may  be  written  in  the  following  form, 

o  =  (hc  +  V  <T,  -  Tref>  <a> 


The  resultant  coefficient  used  in  section  VI  as  h  is  defined  as, 


h  =  h  +  h  . 
c  r 


(b) 


The  formulation  used  for  obtaining  the  coefficients  h£  and  hf  respectively 
are  given  below 


where, 


h  D  i- 

-5 - =  0.  372  Gr  4 

kref 


D  =  diameter  of  the  cylinder, 


=  Grashof  number 
r 


(c) 


krej  =  thermal  conductivity  of  the  gas  surrounding  the  outer  wall 


Equation  (c)  is  now  written  as, 


h  =  k 


ref 


(0. 372)  G~ 


D 


or, 


h  =  k  . 
c  ref 


(0.  372) 
D 


fr,  -  Tref  >  p3  ^ 

T  v2 
1 


(d) 
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where  the  Grashof  number  is  defined  as 


Gr  = 


erefD3 


P  ,=  -  ,  g  =  32.2  — 
T  sec 


and  v  is  the  kinematic  viscosity  of  the  surrounding  gas.  The  term  6 
is  defined  as  the  temperature  difference  between  the  two  surfaces 
exchanging  heat. 

A  further  simplification  of  the  coefficient  h£  is, 


h  =  k  . 
c  ref 


(0.  375) 


(T  -  T  -) 
g  i  ref 


T  v 

i 


The  radient  heat  transfer  coefficient  ,fhr"  used  in  the  problem  is  given 
by  the  formula, 


hr  = 


I  +  _^I 

CI  Att 


f  T:  +  t  Tref  +  T,  Tref  +  Tref  1 
1  108 


where, 

C  =  the  exchange  coefficient  =  108  a  € 

a  =  Stephan  Boltzmann  constant 

e  -  emmisivity 

Aj  =  cylindrical  area 

All  =  surrounding  room  area 

Ai 

Since  the  ratio  of  Aj  to  Ajj  is  small,  the  term  — 

An 

neglected. 


1 
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Now  equation  (f)  can  be  written  as  the  following 


hr 


173 


'  m  3  m2  rp  m  m2  m  3 

i  i  ref  +  i  ref  +  ref 


(g) 


10 


8 


Formulae  d  and  f  added  together  gives  us  the  coefficient  h  which  was 
previously  defined  as 

h  =  h  +  h  .  can  now  be  written  as 
c  r 


h  = 


(0.  372)  (kref) 
1 

D  4 


g  <T  -  Tref)  |  < 


T  v2 
1 


+  0.  173 


T3  +  Ta  Tref  +  T  Traef  +  T^ef 


(h) 


10 


8 


The  standard  formula  used  for  the  over  all  thermal  conductivity  of  the  solid 
layers  of  the  cylinder  (see  figure  7)  is 


K 


m. 


m  \ 

,  / 

D 

D  1 

[_ i  I  + 

~  (in 

3  i  4  1 

—  +  In  — - 

+ 

Id  ) 

K  1 

D 

D  i 

y 

m2\ 

a 

3  1 

' 

(i) 


where  the  thermal  conductivities  are  given  for  the  following  materials 

K  -  dense  partially  stabilized  zirconia 
m 

i 

K  -  cubic  zirconia 
ma 

Km  -  fiberfrax 

11 3 

K  -  steel 
m 

4 
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From  the  shear  relationship,  one  can  arrive  at  the  following  equation  for 


qgw=  0. 176  V  4  P  Cp  „ 


4  lim  fdO  q7~1 
64  o|  d6  J 


where  6  represents  the  radial  difference  (R  -  r)  in  the  fluid  portion  of  the 

1  rR 

cylinder.  The  average  temperature  0  =  -  \  9  2  tt  r  d  r.  (1) 

m  ttR2  ^o 


With  the  help  of  aknown  temperature  distribution  ( 9  ),  h  can  now  be 

D 

computed  anywhere  downstream.  Putting  into  the  formula  for  h  the 

D 

analytical  expression  for  qr  and  #m ,  we  obtain. 
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STABILIZED  ZIRCONIA 


FIRST  EIGENFUNCTION  $.  FOR  d  »0.7 


z 
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SECOND  EIGENFUNCTION  $2  FOR  *  «0.7 


z 
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THIRD  EIGENFUNCTION  $3  FOR*£  «0.7 


0.5  1.0 


FIGURE  4 
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FOURTH  EIGENFUNCTION  $4  FOHaL  »0.7 


smm 


0.5  1.0 


FIGURE  5 


37 


DIMENSIONAL  CROSS-SECTION  OF  THE 
STAGNATION  CHAMBER 


SECTION  A-A 


FIGURE  7 
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EXTERNAL  WALL  TEMPERATURE  DISTRIBUTION 


EXTERNAL  WALL  TEMPERATURE 
DISTRIBUTION  VERSUS  AXIAL  DISTANCE 


AXIAL  DISTANCE- FT. 
FIGURE  9 
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INNER  WALL  TEMPERATURE  DISTRIBUTION  ~#R  X  10- 


INNER  WALL  TEMPERATURE  DISTRIBUTION 

VERSUS  AXIAL  DISTANCE 


FIGURE  10 
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CENTERLINE  TEMPERATURE  DISTRIBUTION 


FIGURE  II 


NON-DIMENSIONAL  INNER  RADIUS- 


RADIAL  TEMPERATURE  PROFILE 


TEMPERATURE  OF  THE  FLUID  -  °RxlO*2 


FIGURE  12 
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RADIAL  TEMPERATURE  PROFILE 
FOR  SPECIFIC  Tw  AND  X 


X*  3.0  FT. 
Tw  =  3565  °R 


TEMPERATURE  OF  THE  FLUID- °R  X  10-2 


FIGURE  13 


NON-DIMENSIONAL  INNER  RADIUS- 


RADIAL  TEMPERATURE  PROFILE 
FOR  SPECIFIC  Tw  AND  X 


X  =  5.0  FT. 
Tw  =3470  °R 
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TEMPERATURE  OF  THE  FLUID- #R  X  10-2 


FIGURE  14 
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NON-DIMENSIONAL  INNER  RADIUS- 


RADIAL  TEMPERATURE  PROFILE 
FOR  SPECIFIC  Tw  AND  X 

XS9.0FT. 


FIGURE  15 
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